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Abstract: By relating the functional averages of a generic scalar operator in simulations 
with Open (O) and Periodic (P) boundary conditions (BCs) respectively for SU{3) lattice 
gauge theory, we show that the scalar glueball mass and the glueball to vacuum matrix 
element can be extracted very efficiently from the former. Numerical results are compared 
with those extracted from the two point function of the time slice energy density (both 
PBC and OBC). The scaling properties of the mass and the matrix element are studied 
with the help of Wilson (gradient) flow. 
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1 Motivation 

In confronting experimental data, lattice Quantum Chromodynamics has achieved remark¬ 
able progress over the years. Nevertheless certain problems emerge as the continuum limit 
is approached, a major difficulty being the spanning of gauge configurations over different 
topological sectors when periodic boundary condition (PBC) is used in the temporal direc¬ 
tion. To overcome this problem, open boundary condition (OBC) in the temporal direction 
has been proposed recently [1-3]. In order to avoid undesirable effects in the spectrum of 
the Hamiltonian, boundary conditions are retained to be periodic for the three-dimensional 
space, which ensures that the transfer matrix is unaltered. Earlier, the advantage of OBC 
over PBC has been recognized and profitably utilized in Density Matrix Renormalization 
Group calculations [4] applied to condensed matter systems. Later, some undesirable fea¬ 
tures of OBC for systems that do not possess an energy gap have been recognized and 
investigated in detail [5]. However, the systems under our consideration, namely, pure 
Yang-Mills theory and QCD fortuitously possess mass gaps. Moreover, OBC yields some 
unexpected extra dividends as we demonstrate in this work. For example, by studying 
the boundary artifacts in the vacuum expectation value of a one point function, one can 
extract the mass and operator matrix elements which are usually extracted from a two 
point function. We illustrate this idea in the context of the calculation of scalar glueball 
mass and glueball to vacuum matrix element in SU(3) Lattice Yang-Mills theory. 

For extracting the mass and the matrix element in lattice Yang-Mills theory, smoothing 
of gauge fields is essential. The Wilson (gradient) flow [6-8] provides a very convenient 
tool for smoothing, with a rigorous mathematical underpinning. Unlike the conventional 
smearing techniques, the Wilson flow provides a common reference scale. Thus by choosing 
a particular flow time, one can study the scaling properties of observables extracted from 
lattice calculations employing different lattice spacings. It is interesting to perform such 
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scaling studies for the glueball mass and the glueball to vacuum matrix element. Another 
noise reduction technique was recently investigated [9] in the extraction of glueball masses. 

2 Relation between correlation fnnctions in OBC and PBC 

We start from the standard Wilson action for SU(3) lattice gauge theory on a x T lattice 
with periodic boundary conditions in all directions 

5'pbc = ^tr[l-Re (2.1) 

^ X fi<U 

where U^y{x) denotes the product of the link variables around a plaquette P in the ~ p 
plane whose lower left hand corner is at x and the sum is over all oriented plaquettes on 
the lattice. 

Using transfer matrix arguments [2, 10, 11], keeping in mind that with open boundary 
condition in the temporal direction, there are no temporal links connecting the time slice 
xo = T — 1 to the time slice xq = 0, one arrives at the action for SU(3) lattice gauge theory 
with open boundary condition 

5'open = ^ [1 - Re U ^, u { x )] (2.2) 

^ X ll<V 

where w{P) is equal to 1 except for the spatial plaquette at time xq = 0 and T — 1 which 
have weight | and here the sum runs over the plaquettes having their corners within the 
time interval [0,r — 1]. 

Thus we find that S'pbc = *S'open + AS where 

AS = ^ I ^ ^ (^tr [1 - Re Uij{y.,T - 1)] + tr [1 - Re [/ij(x, 0)]^ 

^ X I i<j 

+ ^tr [1 - Re 17j4(x, T - 1)1 |. (2.3) 

Note that e~‘^^ is the transfer matrix element [10] between the time slices T — 1 and 
0. Denoting the general transfer matrix element between time slices xq and xq + 1 by 
^-PHm{xo)^ we have 

Hm{xo) = 1 ^ (R “ R® Uij{x,xo)] + tr [1 - Re Uij{x,xo + 1)] ) 

^ X I i<j 

+ ^tr [1 - Re t/i4(x,xo)] |. (2.4) 

An observable which we have extensively studied in our previous work [14] is the time 
slice energy density E{xo) which was used as an interpolating operator to calculate the 
scalar glueball mass. We adopt a symmetric definition of the time slice energy density as 

R (x, xo)] 

^ X I i<j 
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+ ^ X] (ti' [1 “ Ui4{x, xo)] + tr [1 - Re Ui4{x, xq - 1)] ^ | (2.5) 

where C/^;y(x, xq) denotes the oriented plaquette in the fi — i' plane with (x, xq) at its lower- 
left corner. The time slice xq is arbitrary for periodic boundary condition but it is restricted 
within the bulk when open boundary condition is imposed in the temporal direction. For 
the latter case, the definitions of Fl(xo) with xq lying on the boundaries are given by 


B(xo = 0) = 


I i<j 


+ I^tr[l - Re C/j4(x, 0)] I 


( 2 . 6 ) 


and E{xo = T - 1) = ^ tr [1 - Re C/jj(x, T - 1)] 

^ X I i<j 


+ I^tr[l - Re C/j 4 (x,T - 2)] |. 


(2.7) 


To find a relation between (0 (xo))open (^(^o))pbc ''^here O is a generic scalar 
operator we start from 


(O) 


OPEN 


J VU' O{xo) e-^oPBN 

J' VU' 


( 2 . 8 ) 


where the measure VU' excludes the measures for the links connecting the boundary time 
slices (xo = T — 1 and xq = 0). However, as the integrands in both numerator and 
denominator do not depend on these links, we can include their measures without altering 
the result. This leads us to replace VU' by VU which is the measure in case of PBC. Thus 


(0(xo))open ~ 


J VU O{xo) e-'^opEN 
J' VU e“‘S'oPBN 
J VU O{xo) e-^PBc+AS 
j VU e-'S'pBc+AS 


(2.9) 


where the exponents of the integrands in both numerator and denominator on the right 
hand side consist of fields on a periodic lattice. Thus 


(0(xo))open ~ 
= (Cl(a^o))pBC 


J VU O{xo) e-^PBc+A5 1 f VU e-^PBc 
f VU e-'5pBc+A5 ! j VU e-'^PBc 
<0(xo) 


( 2 . 10 ) 


/PBC 
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connected 


( 2 . 11 ) 


= (C>(xo))pBc + ^ (o(xo) 


PBC 


where r = (e^'^)pBC = ^^)pbc- 

As e^^^rn(xo) jg ^^gQ scalai operator, from eq. 2.11 we have 


(C>(xo))open ~ (C’(a:o))pBC + e coshm(^^ - 1 - xq). (2.12) 

where m is the scalar glueball mass. 

In comparison, the two point function for the time slice energy density in the case of 
PBC behaves as 

{E{xo)E{xo = 0))p^^ ^ Co + 2(71 6“”^'^/^ coshm(^^ - xo^. (2.13) 


where 


Cl 


l(0|g(0)|g)l' 

2m 


2m 


(2.14) 


Thus we find that one can extract certain two-point correlators by analyzing the data 
for the functional average of a scalar operator computed with open boundary (in the 
temporal direction) in the region of xq where it differs, due to the breaking of translational 
invariance, from the same computed with periodic boundary. Same technique can, in 
principle, be used to compute any n-point correlator in the scalar channel. In the case of 
lattice QCD with OBC, same technique can also be used to extract the mass of the lowest 
two-pion state. A recent simulation with 2-|-l flavors [16] however encountered large scaling 
violations which unfortunately made such an extraction not possible. 

We note that while the extraction of the glueball mass from the one point function is 
as straightforward as from the two point function, the extraction of the glueball to vacuum 
matrix element from the former is not as straightforward as from the latter. However 
for the operators {E and Elm) used in this work, this becomes possible in the region of 
very small lattice spacing. See the appendix for the discussion. This is supported by our 
numerical results. 


3 Determination of the mass and the matrix element 


In this section, we discuss the methods for extraction of the glueball mass and the glueball 
to vacuum matrix element from the two-point (PBC and OBC) and the one point (OBC) 
correlation functions. 


3.1 Periodic Boundary Condition (PBC) 

In the case of PBC, the mass and the matrix element have to be extracted from the two 
point correlator. Because of the periodicity of the lattice we have 


= G(xo) ~ Co + Cl 


e 


_|_ g-miT-xo) 


{E{xo)E{xo = 0))p3(. 


-mxQ 


(3.1) 





The effective mass is calculated by solving the equation F{m) = 0 


F 


(ri - 1) 


coshm{dt — 1) 


cosh mdt 


+ (1 - r2) 


cosh m{dt + 1) 


cosh mdt 

(3.2) 


where 


n = 


G{xo) ’ 

The coefficient has been extracted as 


— 1) G(a^o +1) ^ rp/^ 

f2 = — -pTt—^dt = T/2 — xq . 


G'(xo) 




where 


= 


aud = _ 

1 Q—ra{xo—l) _j_ g— m(T— xo+l) g—m(xo) g— m(r— xq) 


G(xo) — G(xo + 1) 


g— mxo _j_ g— m,(T—a;o) Q—m{xo+l) g— m(T —xq —1) 

G{xo - 1) - G{xo) 


(3.3) 


(3.4) 


(3.5) 

(3.6) 


3.2 Open Boundary Condition (OBC) 

In the case of OBC, the mass and the matrix element can be extracted from the one point 
correlator as well as from the two point function as discussed in the previous section. In 
the case of two point correlator, because of the lack of translational invariance, we have 


{E{xo + x'o)E{x'o))q^^ = G{xo) Go + 

where Xq is well within the bulk. Thus the effective mass is given by 

, Gjxo - 1 ) - Gjxo) 

G(xo)-G(xo + l)' 

We extract the coefficient as 

c, = i (c« + cf>) 

where 

^( 1 ) ^ G{xo) - Gjxp + 1 ) ^( 2 ) ^ Gjxp - 1 ) - Gjxp) 

1 g—mxo _ g—m(xo+l) 1 g—m(iro —1) _ g—mxo ' 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


4 Numerical results 

SU(3) gauge conhgurations in lattice Yang-Mills theory are generated with open boundary 
condition (denoted by O) at different lattice volumes (by lattice volume we mean total 
number of lattice points) and gauge couplings using the openQCD program [12]. For 
comparison purposes, by implementing periodic boundary condition in temporal direction 
in the openQCD package, we have also generated gauge conhgurations (denoted by P) 
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Lattice 

Volume 

/3 

-^cnfg 

r 

a[fm] 

to jo? 

Oi 

243 X 48 

6.21 

3970 

3 

0.0667(5) 

6.207(15) 

O 2 

323 ^ 04 

6.42 

3028 

4 

0.0500(4) 

11.228(31) 

O 3 

483 X 96 

6.59 

2333 

5 

0.0402(3) 

17.630(53) 

O4 

643 X 128 

6.71 

181 

10 

0.0345(4) 

24.279(227) 

Pi 

243 X 48 

6.21 

3500 

3 

0.0667(5) 

6.197(15) 

P 2 

323 ^ 04 

6.42 

1958 

4 

0.0500(4) 

11.270(38) 

P 3 

483 X 96 

6.59 

295 

5 

0.0402(3) 

18.048(152) 


Table 1. A^cnfg is the number of configurations, r is the trajectory length used in the HMC 
simulation and is the dimensionless reference Wilson flow time. O and P refer to ensembles 

with open and periodic boundary conditions in the temporal direction. 


for several of the same lattice parameters. Simulation details are given in table 1. The 
parameter to is defined in the context of Wilson flow [6-8] which is used to smooth the 
gauge configurations. The implicit equation 

{t\E{T/2))}^^^^=0.3 (4.1) 

with t and T being respectively the Wilson flow time and the temporal extent of the 
lattice, defines a reference flow time to which provides a reference scale to extract physical 
quantities from lattice calculations. The effectiveness of the Wilson flow in the extraction of 
topological susceptibility [13], glueball mass [14] and topological charge density correlator 
[15] has been demonstrated recently. 


0.0026 


0.0024 

& 0.0022 

0.002 

0.0018 


Figure 1. Plot of {E{xo)) versus Xq at flow time t = to at (3 = 6.59 and lattice volume 48^ x 96 
for ensemble O3 (Hlled circle) and ensemble P 3 (filled square). 

In hg. 1 we plot the configuration average of the time slice energy density averaged 
over the spatial volume {E{xo)) versus the time slice xo at the reference flow time t = to 


P = 6.59 
V^= 0.35 fm 


“ O3 

. Po 
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[7] at /3 = 6.59 and lattice volume 48^ x 96 for ensemble O 3 (filled circle) and ensemble P 3 
(filled square). We find that in the boundary region at xq = 0 , (P(xo)}oBC rises above and 
then decays to (i7(xo))pBC- Similar behaviour is found in the boundary region xq = T — 1. 
As explained in sec. 2 , the decay rate is determined by the mass of the lowest excitation 
in the scalar channel, namely the glueball mass. 



Figure 2. Plot of tg(£'(xo, to)) (left) and to(i7m(a;o, to)) (right) versus xo/i/to for the ensembles 
Oi, O 2 , O 3 and O 4 . 

In order to reliably determine the glueball mass from the decay of (P(xo)}oBC and 
(ffm(xo)}oBC first we need to verify their scaling behaviour. Towards this goal, in hg. 2 
(left) we plot tQ(E(xo,to)} versus xo/-v/to for the ensembles Oi, O 2 , O 3 and O 4 . We note 
that, except for the largest lattice spacing, the data exhibit excellent scaling behaviour in 
the tail region from where one can extract the glueball mass. In comparison, in hg. 2 (right), 
the scaling behaviour of tQ{Hm{xo,to)) is shown. The difference between these observables 
decreases as the lattice spacing decreases. The slightly worsened scaling behaviour of the 
latter for the two relatively larger lattice spacings can be readily attributed to the following. 
The expression for Hm follows directly from the lattice action used in this work, namely, the 
unimproved Wilson gauge action. On the other hand, the expression for E uses the clover 
definition of the lattice field tensor. Thus the difference in scaling behaviour exhibited in 
the left and right parts of fig. 2 results from the use of an improved versus an unimproved 
operator, which diminishes as lattice spacing goes to zero. The reason for the non-smooth 
behaviours of both E and Hm in case of the ensemble O 4 is due to the lack of statistics. 

In fig. 3 we show an example of the effective glueball mass determination from the one- 
point function. The effective mass m(0'’”*“) is plotted versus xq for four different Wilson 
flow times for the ensemble O 3 . We find that there is an optimum window of Wilson 
flow time within which the glueball mass can be reliably extracted from the effective mass 
plot. For lower values of Wilson flow time, the smearing is not able to remove unwanted 
contributions completely and the plateau region is too narrow. For higher values of Wilson 
flow time also plateau is too narrow, presumably due to over smearing resulting in the 
overlap of the two glueball interpolation operators in the correlation function. Ideally the 
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Figure 3. The scalar Glueball effective mass as a function of Xq extracted from the 

one-point function of E for four values of Wilson flow times for the ensemble O 3 . 

plateau region should be independent of the flow time chosen, but we observe that an 
overlapping (within the statistical errors) and extended plateau exists only in the region of 
flow time between 0.3 fm and 0.35 fm. In the following we will present the glueball mass 
and glueball to vacuum matrix element extracted for these two values of the flow time. 

In table 2 we present the values of the lowest scalar glueball mass me (= m/a) in MeV 
and the glueball to vacuum matrix element in unit of the Sommer parameter ro extracted 
from correlators at Wilson flow times = 0.3 fm and 0.35 fm. One point refers to 
the one point function of Hm and E and Two point refers to the two point correlator of 
E. Wilson (gradient) flow is known [6-8] to have the consequence that the expectation 
values of local gauge invariant operators constructed from the gauge field at positive flow 
time are ultraviolet finite. Thus we expect the glueball mass and the glueball to vacuum 
matrix element extracted at a fixed flow time at different lattice spacings to exhibit scaling 
(provided lattice artifacts are negligible). The glueball mass should be independent of the 
flow time but the glueball to vacuum matrix element is expected to depend on the flow 
time (energy scale). The results presented in table 2 are consistent with these expectations 
within the statistical errors and the limited range of flow times probed. The investigation 
of the relationship between the extracted glueball to vacuum matrix element at a given 
Wilson flow time and its continuum counterpart involves a detailed numerical study of the 
behaviour of the glueball matrix element of the energy momentum tensor and the trace 
anomaly under Wilson flow (for related theoretical work, see for example. Refs. [17] and 
[18]) and is beyond the scope of the present work. 

In fig. 4, we have presented the variation, with a^, of lowest glueball mass and glueball 
to vacuum matrix element in unit of Sommer parameter ro extracted from one-point (Open) 
and two-point correlators (Open and PBC) of E for different lattice spacings and lattice 


- 8 - 














(ASlAl) 


Lattice 

Correlator 

Mass mo (MeV) 

Coefficient r^C 

\/8t=0.3fm 

\/^=0.35fm 

-v/^=0.3fm 

v/^=0.35fm 

Oi 

One point (Hm) 

1729 (67) 

1649 (56) 

98 (8) 

85 (6) 

One point (E) 

1690 (94) 

1650 (80) 

104 (14) 

95 (11) 

Two point 

1626 (186) 

1501 (96) 

103 (17) 

88 (7) 

Pi 

Two point 

1625 (92) 

1594 (72) 

106 (10) 

100 (7) 

O 2 

One point (Hm) 

1700 (70) 

1645 (59) 

102 (8) 

90 (6) 

One point {E) 

1710 (84) 

1629 (48) 

113 (12) 

97 (5) 

Two point 

1587 (234) 

1458 (115) 

101 (21) 

85 (8) 

P 2 

Two point 

1552 (61) 

1506 (71) 

99 (5) 

91 (6) 

O 3 

One point(Lfm) 

1640 (101) 

1551 (81) 

103 (14) 

88 (9) 

One point {E) 

1625 (85) 

1540 (69) 

108 (12) 

93 (9) 

Two point 

1616 (254) 

1465 (161) 

99 (20) 

84 (12) 

P 3 

Two point 

1467 (181) 

1421 (182) 

94 (16) 

85 (16) 

O 4 

One point {Hm) 

1818 (87) 

1752 (177) 

140 (13) 

125 (27) 

One point {E) 

1783 (141) 

1711 (273) 

144 (24) 

126 (48) 

Two point 

1521 (513) 

1459 (488) 

90 (35) 

83 (33) 


Table 2. The lowest scalar glueball mass (ma) in MeV and the glueball to vacuum matrix element 
in unit of the Sommer parameter tq extracted from correlators at Wilson flow times ^/M — 0.3 fm 
and 0.35 fm. One point and Two point refer to the one point and the two point correlator of E. 

200 


150 


u 

m 100 
o 

50 


0 

0 0.001 0.002 0.003 0.004 0.005 0 0.001 0.002 0.003 0.004 0.005 

(fm^) (fm^) 

Figure 4. Plot of lowest glueball mass in MeV (left) and glueball to vacuum matrix 

element in unit of rg (right) versus extracted from one-point (Open) and two-point correlators 
(Open and PBC) of E for different lattice spacings and lattice volumes at Wilson flow time -s/St = 0.3 
fm. Light blue and light violet shaded regions correspond to the error bands around the fit curves 
to one-point and two-point data respectively. 




volumes at Wilson flow time \/^ = 0.3 fm. Also shown are separate fits to one-point and 
two-point data sets. Light blue and light violet shaded regions correspond to the error 
bands around the fit curves to one-point and two-point data respectively. Corresponding 
data at Wilson flow time \/^ = 0.35 fm are plotted in fig. 5. 
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Figure 5. Plot of lowest glueball mass in MeV (left) and glueball to vacuum matrix 

element in unit of rg (right) versus extracted from one-point (Open) and two-point correlators 
(Open and PBC) of E for different lattice spacings and lattice volumes at Wilson flow time VM = 
0.35 fm. Light blue and light violet shaded regions correspond to the error bands around the fit 
curves to one-point and two-point data respectively. 


Correlator 

Mass me (MeV) 

Coefficient r^C 

\/^=0.3fm 

\/^=0.35fm 

\/^=0.3fm 

\/^=0.35fm 

One point {Hm) 

1701 (44) 

1628 (36) 

100 (5) 

88 (4) 

One point {E) 

1674 (51) 

1610(35) 

108 (7) 

95 (4) 

Two point 

1572 (46) 

1490 (61) 

100 (4) 

90 (3) 



1500 


1000 - 


500 - 



PBC 2-point I— •- 
Open 2-point i— ©- 
Open 1-point i— A- 

Fit to 2-point data - 

Fit to 1-point data - 


Table 3. Fit results for glueball mass and glueball to vacuum matrix element. 


Note that with PBC, glueball observables can be extracted only from the two point 
correlator. With PBC, the signal can be extracted after performing source averaging and 
hence the statistical error can be reduced. However, when one simulates at lower lattice 
spacings it becomes increasingly difficult to generate statistically independent configura¬ 
tions with PBC in the temporal direction and hence the statistical accuracy begins to 
suffer. This problem is overcome with OBC in the temporal directions. However in this 
case translational invariance is lost close to the boundary and hence one can perform source 
averaging only over the spatial volume and a few temporal slices well inside the bulk. This 
effectively increases the statistical error. With two point correlator, the disconnected con¬ 
tribution has to be subtracted numerically and hence this contributes to the increase of 
statistical errors since one is dealing with the subtraction involving two large quantities. 

We have shown that, with OBC, the glueball observables can be extracted alterna¬ 
tively from the one point correlation function as well. The automatic subtraction of the 
disconnected part in the case of one point function leads to smaller statistical error com¬ 
pared to the extraction from two point correlation function with same boundary condition 
as exhibited in table 2 and figures 4 and 5. One should, of course, keep in mind that this 
method for the calculation of n-point correlators involving the energy density is applicable 
only in the scalar channel. 
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It is noticeable that the results (for both mass and matrix element separately) extracted 
from the two point correlation functions with periodic and open boundary conditions are 
very close to each other. However, there is systematically an upward shift in the results 
obtained from one point correlation functions (doable only with open boundary condition), 
although not beyond the statistical errors in most cases. We have treated the data sets 
extracted from two point and one point correlation functions separately in the fitting 
procedure. It is also noteworthy that for the range of lattice spacings explored in this 
work, scaling violations are within the statistical uncertainty of our data. This led us 
to perform just constant fits to the data. While fitting the data (both mass and matrix 
element) for one point correlator, we have excluded f3 = 6.71 as it somewhat deviates from 
the general trend. However, we have checked that its inclusion in fitting procedure does 
not change the results significantly because of large errors on the data at this coupling. 
The fit results are presented in table 3. 

In our previous work on glueball mass extraction [14], a systematic study of the varia¬ 
tion of mass with Wilson flow time was not performed. For each (3, flow time yielding the 
most stable plateau was picked up. However, for the study of glueball to vacuum matrix 
element it is mandatory to choose a common flow time for all the lattice spacings. This 
path has been followed in the present work. This causes the difference in results for glue¬ 
ball mass obtained via the two point correlation functions in this work from that quoted 
in [14]. However, it is gratifying to note that the average of the glueball masses obtained 
with the two flow times is very close to the value quoted in [14] and also agrees well with 
the value given in [19]. Although our results for glueball mass extracted from one point 
correlation function are comparatively higher, they fall between the results of [19] and [20] 
both evaluated using two point correlators. 

Conclusions 

In lattice QCD, OBC in the temporal direction has been proposed to overcome the 
difficulty in the spanning of gauge configurations over different topological sectors. How¬ 
ever, the lack of translational invariance in this case can cause some inconveniences in the 
measurement of observables due to boundary effects when compared to the case of PBC. 
In this work, we have demonstrated that the same boundary artifacts can be exploited to 
yield certain observables with greater efficiency. This is achieved by relating the functional 
average of a generic scalar operator measured in the case of OBC to that with PBC. The 
scalar glueball mass and the glueball to vacuum matrix elements obtained from this ob¬ 
servable in the case of OBC are compared with the values extracted from the measurement 
of two-point function of the time slice energy density in the case of both PBC and OBC. 
The Wilson (gradient) flow is used to exhibit the scaling properties of both the mass and 
the matrix element. 
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A Connection between matrix elements 


In the following, we show that, with a set of approximations, one can understand why 
the glueball to vacuum matrix element of both and E calculated from the one point 
functions agree (within statistical errors) with that calculated from the two point function 
of E. 

Consider the two point correlator appearing in eq. 2.11 for the one point function of 
E{xo). In this correlator, for xq, in E(xo), far away from the boundaries the last term in 
the expression for AS* given in eq. 2.3 can be approximated as 
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5" 


^ J^tr[l-Re [/i 4 (x,r-l)] 
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X L 


2 I ^ Ci 4 (x, 0)] ^ tr [1 - Re 17*4(x, T - 1)] 


+ ^ tr [1 - Re t/*4(x, T - 1)1 + ^ tr [1 - Re [/*4(x, T - 2)] ) 1. 

i i ) 

With this approximation we find, A5 ^ ^77(0) -|- E{T — 1)^. This leads us to write 
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Now in continuum limit E{xo) = aEyMixo) where Eym{xq) = 1^ {G^iu (x) ( 2 :)}. 

X 

This enables us to write 


^ 'e(O) -I- E(T-1) 


— 1 + a—i^EYM{^) + Eym{T — + O . (A.2) 


So eq. A.l can be approximated as 
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As evident from fig. 2, Hmixo) can be approximated by E{xq) in the limit of very 
small lattice spacing. Therefore, in this limit one point function of takes the form 


-ected 



{Hm{xo)) 


OPEN ~ 



(A.4) 


Discussions subsequent to eq. A.l upto eq. A.3 follow thereafter. 
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